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Recent results relevant for the B-physics phenomenology, obtained from lattice QCD simulations by the APE 
Collaboration, are reviewed. This includes the B° — B° mixing amplitude, B — > ty semileptonic decay and the 
relative width difference of B® mesons, (Ar/r)s s . 



The main theoretical obstacles in determining 
the amount of CP-violation that comes from the 
Standard Model (SM) are related to the uncer- 
tainties in computation of various hadronic quan- 
tities. In this talk, I focus on several such quan- 
tities/processes involving heavy-light mesons, for 
which APE group provided new lattice results. 
Technical details about the simulations are given 
in the references which will be quoted with each 
quantity I discuss in this paper. Here I only stress 
that all our results are obtained at (3 = 6.2, in the 
quenched approximation, and by using the (fully 
propagating) 0{a) improved, Wilson fermions. 

1. Decay Constants (fg and others) [ |^] 

One of the essential hadronic quantities enter- 
ing the — mixing amplitude is the B-meson 
decay constant fg (q — d,s), defined as 

(0\h^ 5 q\B q {p)) = ipjg t . (1) 

The central results from our two simulations (see 
Tab. |l|) are obtained by: (i) linearly extrapolat- 
ing (interpolating) /// ?v /m^7 in l/m Hq , to the B 
(D) meson mass; (ii) including the KLM factor 
in a way discussed in [ |J| ; (iii) converting to the 
physical units by using a~ 1 {m,K*) — 2.7(1) GeV. 
To estimate the systematic errors, we combine in 
quadrature the following differences between our 
central values and the ones obtained when: (a) 
extrapolating in l/m,H q quadratically; (b) ommit- 
ing the KLM factor, (c) using the ratio fg/fn 
to extrapolate to mg. Since the extrapolation 
from the directly accessed heavy meson masses 
(2 GeV < mg < 3 GeV) to mg is rather long, 
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the source (a) largely dominates the systematics. 
This error has not been included in the results [ , 
where only 3 heavy quarks were considered. Note, 
however, that this error completely cancels in the 
SU(3) breaking ratio, fg s /fg d . For comparison 
with results of other lattice groups, see [ ||. In [ 



Table 1 

Results for pseudoscalar decay constants in MeV. First er- 
rors are statistical and the second are systematic. 





Ref. [ 




I 


Rcf. [ 2 




f B 

fB s 1 JB 

Jd 
fo s 1 Id 


173(13)+f 
1.14(2)(1) 
216(H)(5) 
1.11(1)(1) 


175(22)^ 
1.17(4)(1) 
207(ll)±o 
1.13(3)(1) 



|3| we have also computed the vector meson decay 
constants, which I now update. In addition we 
compute the coupling of the vector meson to the 
tensor current, i.e. 

(0|6vz|B*(p,A)) =ie^m Bs f B; 
(0\ba^q\B* q (p,X)) =i{e^p v -p^)fl«{ii) , 

where fi is the scale at which the tensor density is 
non-perturbatively renormalized (in the RI-MOM 
scheme). These decay constants are particularly 
important in testing the validity of the factoriza- 
tion in non-leptonic decays of heavy-light mesons. 
Our new results are 

fg* = 199(14)+f MeV ;f D , - 258(14)(6) MeV; 
/u;//u. = 1-14(2) (1) ;fl.{m h )/fg. = 0.88(2)(2) ; 
fD*Jfn. = 1-10(2) (2 GcV)/f D * = 0.90(2) . 
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2. B°-B° Mixing and (AF/F) Bs [ g g§ 

To access any bare continuum Ai? = 2 opera- 
tor from the lattice, by using Wilson fermions, one 
first has to subtract the effect of mixing with other 
dimension-six 4-fermion operators which is due to 
the explicitely broken chiral symmetry in the Wil- 
son quark action. A bare (continuum) operator 
should then be appropriately renormalized. The 
whole procedure can be shortly written as 
(B Q q \Q^)\B q ) = 



(5o;M) 




where k run over the basis of parity conserv- 
ing operators (Qi), Ai(go) and 2ij(0o>M) are the 
subtraction and renormalization (RC) constants, 
respectively. A technique to compute the con- 
stants Ai and Zij non-perturbatively, in the RI- 
MOM renormalization scheme, has been developed 
in ref. [ |5|. We work in Landau gauge, apply the 
technique [ ||] at three different values of the scale 
fx, and verify that the scale dependence of Zij(fi), 
for the operators discussed below, is indeed well 
described by the perturbative NLO anomalous di- 
mension matrix [ |(| . This allows us to express our 
matrix elements in the renormalization group in- 
variant (RGI) form. 

B q —B q Mixing: The needed parameter, B Bq 
(q = d, s), is defined as 

(B° q \Q L (ri\B° q ) = \mlj% q B B », (2) 

where Q L = 6*7^(1 - 75)<? 1 b 3 !^ ~ lb)q J , and i,j 
are the color indices. From the definitions (0) and 
((sj) , it is clear that the _B-parameter can be di- 
rectly accessed if (for each heavy light-meson H q ) 
we compute the ratio of correlation functions 



BhM ( 3 ) 
The last limit is valid when the operator Ql and 
the pseudoscalar sources P are sufficiently sepa- 
rated on the temporal axis. To reach the physi- 
cally relevant B Bq from the extracted Bjj q (which 

2 The procedure sketched above can be highly simplified if 
one uses the Ward identity to relate the parity conserving 
to the parity violating operators (for which = 0) 
This idea is yet to be implemented in practice. 



should scale with heavy meson mass as a constant) , 
we make the linear \jm B fit and extrapolate to 
m B . Our final results are 



B Bd = 1.38(11) 



+ .00 
-.09 ) 



Br 



1.35(5) 



+ .00 
-.08 ' 



B B ,/B Bd =0.98(5) 



(4) 



where we converted (to NLO) the directly com- 
puted B™- MOM {[i) to the RGI form B Bq ■ We also 

obtained, B D = 1.24(4)t ™, which may be useful 
in the non-SM phenomenology. Since we clearly 
see 1 corrections from our data, one can try 

to constrain the extrapolation by using the static 
result for B B [M. Such an exercise leads to a ~ 5% 
lower value of B B , which is well within our error 
bars (a similar conclusion is reached in [ Q). 

After combining the above results with the ones 
for f B , also computed in [ |j , we get 




f Bd \/B Bd = 206(28) (7) MeV; 
I~B~r 

1.16(7) 



(5) 



In the last result, most of the systematic uncer- 
tainties cancel in the ratio. To exemplify the 
phcnomcnological benefit of this result, I combine 
our value for £ with the updated experimental 
value for the frequency of B® mesons oscilations 

(Amf p,) = 0.486(15) ps^ 1 [ §), to get 



Am. 



\Vts\ 2 £ 2 ( m B s 



. , v , Am d 

\VtdV \ m B d 

16.2 ±2.1 ±3.4 ps -1 , 



(cxp.) 



(6) 



where |V ts | 2 /|V t d| 2 = 24. 4±5.0 is assumed. Exper- 
imental lower bound is Ami° xp '' 1 > 14.9 ps -1 
(AF/r) B : In the framework of the heavy quark 
expansion, the leading contribution in the expres- 
sion for the width difference of B® mesons, comes 
from AB — 2 operator, Q s = b l {l - 75)5*62(1 - 
75)s J , the matrix element of which is parameter- 
ized as 

5 / m%J Bs 
3 \m b (fi) + m s (/i). 



(B° s \Qs^)\B° s ) 



where i?s(/i) is the wanted bag parameter. An im- 
portant observation made in ref. [ [l^] is that if we 
only replace Ol/8 — > — Os/5 in eq. (||), we see a 
very large dependence on the inverse heavy meson 
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mass l/jriH,- On the contrary, if in denominator 
of eq. (g) we also replace the axial current by the 
pseudoscalar density, Aq — > P, the l/rriH s depen- 
dence becomes much weaker and the extrapolation 
to m B is more under control (which is why our cen- 
tral results are those obtained using the latter pro- 
cedure). Obviously, the large 1/mn dependence 
comes from the ratio of the heavy meson/heavy 
quark mass, [m^/fmg + m s )] 2 . It will be inter- 
esting to see whether the inconsistency of the two 
procedures disappears with the simulations per- 
formed closer to m Bs . In this calculation, we also 
needed to compute the matrix element of the oper- 
ator Qs which is the one obtained by reversing the 
color inidices in Qs, and which mix with Qs in the 
continuum. Once we extract the values for -Bs(/i) 
and B s {p), we converted RI-MOM^ MS since the 
Wilson coefficients, G{z) and Gs{z) (z — ml/ml), 
were computed in the MS scheme [ |ll| . After a lin- 
ear extrapolation to 1/mg,, we get 

B™(rn b ) = 0.86(2)(3); = 1.25(3)1;°* . 

For the physical prediction of (AT/T) B , one needs 
the ratio of the matrix elements (Q) and (Q). We 
obtain 



ft Mb (^ = m b ) = -0.93(3)(l) 



(8) 



which is in a good agreement with results obtained 
by using the effective theories [ |§|, [l2]| . We pro- 
posed in [ [l(]] a safer way to predict 



K T Bs Am Bd 



m Bd 



(oxp. 



G(z)-Gs(z)K(m b ) + ~6 1/m ) ? (9) 



where K is the known constant and 5i/ m encodes 
l/m b corrections (which are estimated by using 
the vacuum saturation approximation). The ad- 
vantage of using this formula is that it contains 
experimentally well determined quantities, and £ 
and lZ(m b ) ratios, in which (again) most of the sys- 
tematic uncertainties cancel. Finally, we obtained 



AT 



= [(0.5 ±0.1)- (13.8 ±2.8)K{m b ) 

+(i5.7±2.8)5 1/m ] -nr 2 

= (4.7 ± 1.5 ± 1.6) • 10 -2 , (10) 



where we show how the explicit cancellation occurs 
between the leading lZ(m b ) and the l/m b correct- 
ing terms. Therefore, to improve the above result 
it is necessary to gain a better control over the 
dimension-seven operators which appear in 5\/ m . 

3. | Kb | from B(B -> irtv) [ |]| 

Compared to "Lattice 99" [ [jj , we now have the 
final results for the D decay modes (see [ and 
also for the most challenging one, B — > tt. The 
relevant form factors, F + / (q 2 ) (q = p B —p), 



(*(p)\M\H(p H )) 



^Mq 2 ) 



Ph +P 



H 



<r 



F+{q 2 ) , (11) 



are extracted for 3 different light (q) and 4 heavy 
(Q) quark masses and for 13 inequivalent kinemat- 
ical configurations (J>H,q)- The mass extrapola- 
tions of form factors are known to be trickier be- 
cause of the interplay between ra B and q 2 depen- 
dences. A parameterization for the q 2 -dependence, 
which contains most of the theoretically available 
constraints, has been proposed in [ Q : 

c H (l - a H ) 



W) = 



F (q 2 ) = 



(1 - q 2 ) (1 - a H q 2 ) 

cg(l - ttjf) 
(1 



(12) 



where q 2 



9 I 

q jm 



The parameters </> G 
{cn^/rriH, (1 — an)mH, (0 — l)m_f/} should scale 
as a constant (plus corrections oc l/mn)- To reach 
B — > 7r we first fit the form factors to the param- 
eterization ( |l2[ ) for each combination of the heavy 
and light quarks and then adopted the following 
two methods: 

Method I We smoothly extrapolate to the final 
pion (kaon) state, for every heavy quark, and then 
use the scaling laws for all three parameters to ex- 
trapolate to B, namely 

4> = a H 

Method II Following the proposal of UKQCD 
group [ [l^] , we first extrapolate to the final pion 
at fixed 

Ml-q 2 



Q2 

2 
H 



in 



(13) 



v ■ p 



Ml 



2M 



H,, 



(14) 
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taoie z 

APE lattice results for the B — > 


7r form factors [|l3) are 


compared to the predictions obtained by using 


LCSR [|T|. 




Method I 


Method II 


LCSR [ 


16| 


CB 


0.42(13)(4) 


0.51(8)(1) 


0.41(12) 




0.40(15)(9) 


0.45(17)l;°j 


°- 32 -.07 


(3b 


1.22(14)t^ 


1.20(13)1$ 






F(0) 


0.26(5) (4) 


0.28(6)(5) 


0.28(5) 


F B ^ K {0)/F B ^(0) 


1.21(9)1;°° 


1.19(11)1;°? 


1.28(11)" 


f.18 
-.10 


F (m%) 


1.3(6)1;° 


1.5(5)1;° 






gB'B-R 


20 ±7 


24 ±6 


22 ±7 


|K6|- 2 r(B° -» n-l+v) [ps 


- 1 ] 6.3 ±2.4 ±1.6 


8.5 ±3.8 ±2.2 


7.3 ±2 


.5 



and then extrapolate to B by using the HQET 
scaling laws. The g 2 -behavior is then fit by using 
eq. @. 

The results obtained by using both methods are 
shown in Tab. [^, where we also make a comparison 
with the light cone QCD sum rules (LCSR) predic- 
tions [ |l6| . Besides a very good consistency of the 
results of the two (different) methods, a pleasant 
feature of this analysis is also the agreement with 
lattice results of [ [l5| as well as with the LCSR re- 
sults [ [l^] . For a further comparison with presently 
available lattice results, see [ [|. 

The central values in Tab. || are obtained 
through a quadratic extrapolation to tub, which 
provides a better consistency with the Callan- 
Treiman relation, F (m 2 B ) = fs/fnU Tnc 
complete account of the systematic uncertain- 
ties is detailed in [ |l3). I emphasize that 
F B ^ T (0) in Tab. ||, as obtained after extrapolat- 
ing F H ^™ (0)m 3 fl 2 to wig, is indistinguishible from 
the one obtained by combining separately extrap- 
olated Cb and as into ^^^(0) = Cb(1 — qb). 

To read the content of the parameterization ( p"2| ) 
I recall that the parameter cb measures the residue 
of F+(q 2 ) at m 2 B ,, thus allowing determination of 
the qb* Btt • By using the standard definitions, one 
has gB'B-K = 2ms* cs / '/b* and the results are 
also given in Tab. ^. The parameter as indi- 
cates that besides the first pole (q 2 = m 2 B ,), all 



3 Research on verification of this relation on the lattice has 
been discussed at this conference [ h?j] . 



the other singularities contributing to F + {q 2 ) can 
be mimicked by an effective pole corresponding to 
m\- ~ 8 ± 1 GeV, whereas the contributions to 
the scalar form factor Fo(q 2 ) are represented by 
an effective pole with the mass m + ~ 6 ± 3 GeV. 

Finally by comparing our results for the de- 
cay width, with the experimental branching ratio 



B{B° 



-£D), we obtain \V ub \ = (4.1 ±1.1)10" 
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